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AN IMPROVED A PRIORI ERROR ANALYSIS OF NITSCHE’S 
METHOD FOR ROBIN BOUNDARY CONDITIONS 

NORA LUTHEN, MIKA JUNTUNEN, AND ROLF STENBERG 


Abstract. In a previous paper [6] we have extended Nitsche’s method [8] for 
the Poisson equation with general Robin boundary conditions. The analysis 
required that the solution is in H'^, with s > 3/2. Here we give an improved 
error analysis using a technique proposed by Gudi [-j]. 


1. The method and its consistency 


In the article [6] a Nitsche-type method is introduced and analyzed for the 
following model Poisson problem with general Robin boundary conditions: Find 
u G H^{n) such that 


(1.1) —Au = f in n, 

du 1 

(1.2) -^ = -[uo-u)+g on r, 


where fl C N = 2,3, is a bounded domain with polygonal or polyhedral 
boundary F, / G L^(fl), uq G il^/^(r), g G L^(r), and e G K, 0 < e < oo. 
The limiting values of the parameter e give the Dirichlet and Neumann problems, 
respectively. 

The error analysis presented was not entirely satisfactory. It assumed that the 
solution is in 77^(17) with s > 3/2, which is the same condition that traditionally has 
been needed for discontinuous Galerkin methods [4]. For discontinuous Galerkin 
methods Gudi introduced a technique using a posteriori error analysis by which 
this assumption could be avoided [-5]. 

The purpose of this paper is to use these arguments to improve the analysis 
of the Nitsche method for the above Robin problem. Below we start by recalling 
the method of [6]. We first recall the derivation of the method in a way that 
emphasizes the use of the residual, which will be crucial for the error analysis. The 
same notation as in [6] will be used. The finite element partitioning into simplexes 
is denoted by 7/. This induces a mesh, denoted by Gh, on the boundary F. By 
K G Th we denote an element of the mesh and by E we denote an edge or a face in 
Gfi- By hji we denote the diameter of the element A G 7/, and by px the radius 
of the biggest ball contained in K. The mesh is assumed to be regular, i.e. it holds 


(1.3) 


hx 

sup — = K < oo. 
Ken Pk 


By we denote the diameter oi E G Gh- The finite element subspace is denoted 
by 

Vh := {v G : v\x G iPp(A) VA G % }, 

where VpiK) is the space of polynomials of degree p. 
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The Nitsche method is obtained as follows. Multiplying the differential equation 
(1.1) with a testfunction w € and integrating by parts we have 


(1.4) 

Defining the residual 

(1.5) 

the boundary condition is 

( 1 . 6 ) 

Hence it holds 

(1-7) E 


, du 


{Vuyw)^-{ — ,w)^-{f,w)^ = 0. 


Rr{v) = e(|^ - g)+v-Uo, 


Rt{u) = 0. 


1 


EeQ, 


e + 'yhE 


{Rriu),w)j^ = 0 


and 

( 1 . 8 ) 


-E 

E^Qh 


ihE , , . dw. _ 

-{Rr{u), 0 . 


e + y/iE 


Adding (1.4), (1.7) and (1.8) shows that the exact solution satisfies 
, du 
'dn 

(1.9) 




EGGh 

- E 


e + y/iB 


jhE /p / ^ _ 

[Rriu), 0. 


Eeg^ 


e + 'yhE 


Substituting the expression (1.5) for the boundary condition and rearranging the 
terms, we see that the exact solution satisfies 


( 1 . 10 ) 

where 


Bh{u,w) - Th{w) = 0 \/wGVh 


EeQh I ' ^ 


1 , . EyhE /dv dw, 

^ - ^ + ^hE 


e + 'yhE 


and 


(1.12) Mw)= {f,w)^+ E I 

EeGh I 


1 , , 'fhE / dw , 


e + 'yhE 


dn ' 


e + 'fhE 

The above derivation shows the consistency of the 
Nitsche Method [6]. Find Uh G 14 such that 
(1.13) Bh{uh,w) = Fhiw) yweVh. 


e , . ejhE / dw. 


e + 'fhE ' ’ dn'^ \ 
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2. The new a priori error estimate 


The estimate will be given in the mesh and problem dependent norm 

(2.1) lkllI:=l|VT||g,o+ E 

We recall the following discrete trace inequality which is easily proved by scaling 
arguments. 


Lemma 2.1. There is a positive constant Cj such that 

2 

< C/|lVz;||2_o ^veVh. 

0,E 

For the formulation we have the following stability result, cf. [ 6 ]. Here and in 
what follows C denotes a generic positive constant independent of both the mesh 
parameter h and the parameter e. 

Lemma 2.2. Suppose that 0 < 7 < I/C 7 . Then there exists a positive constant C 
such that 

(2.3) Bh{.v,v)>C\\v\\l VvGVh. 

By fh S Vh and gu, UQ^h G 14|r we denote the interpolants to the data. For 
E G Qh we denote by K{E) G Th the element with E as edge/face. In [ 6 ] we proved 
the following bound. 


( 2 . 2 ) 




dv 

dn 


Lemma 2.3. For an arbitrary v G Vh and E G Gt it holds 


l 1/2 

(2.4) ^^|li?r(i;)||o.E < C(|| V(u - u)||o,7c(i5) + /ik||/ - Mlo.ic(E) 


6 + hi 


1/2 


+ 7 L “ ^lloT + l|e(g - 9h) + Uo - Mo.ftllo.u). 

(e + /iE)u2 ’ ’ / 

We introduce the oscillation terms 

1^2 

(2.5) osc(/) = [Y. - MIk) , 

KeTh 

( 2 . 6 ) osc{e,uo,g) = ( n 2 lk(g " g^) + - mo.IiIIo.e) ■ 


Lemma 2.3 then gives 
Lemma 2.4. For v GVk it holds 

(2.7) ( Y ( Ih w ll-^r(i^)||o,E) ^ < C'{||M-^;|4+osc(/) + osc(e,uo,5)}• 

Et^. ’ 

We can now prove our new error estimate. 


Theorem 2.1. Suppose that 0 < 7 < IjCj. Then there exist a positive constant C 
such that 

( 2 . 8 ) \\u - Uh\\h < C{ inf ||u - u|4 + osc(/) + osc(e,ito,5)}. 

veVh 
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Proof. We will divide the proof in 6 steps. 

1. Treating the consistency by Gudi’s method. 

Let V G Vh he arbitrary. From the stability we have 

(2.9) C\\v - UhWl < Bh{v - Uh,v - Uh). 

Next, we denote w = v — Uh and use (1.13) 

(2.10) Bh(v - Uh,v - Uh) = Bh{v - Uh,w) = Bh{v,w) - Bh{uh,w) 

= Bh{v,w) - Fh(w). 


Reversing the arguments leading from (1.9) to (1.10) we see that 


dv 


Bh{v,w) - Fh{w) = { — {f,w) 


( 2 . 11 ) 


+ E 

EeOh 

- E 

EeOh 


1 


e + 7^E 


{Rt{v),w)^ 


ihs /„ , , dw. 


Substituting the boundary condition (1.2) into (1.4) we get 


( 2 . 12 ) 


{Vu,Vw)^ - {-{uo - u) +g,w).^ - {f,w)^ = 0. 


Subtracting this from the right hand side of (2.11) yields 

, dv 1 
' dn e 


Bh{v,w) - Fh{w) = {V{v - u),Vw)^---{uo - u) - g,w)^ 

1 


+ E 


EeGh 


e + -fhE 


{Rr{v),w) 


E 


(2.13) 


- E 


ihE /p , s dw. 


e + 'yhs 


EgQh 

= Ri + R2 H" R'^ H“ Ra- 


dn' ^ 


Next we estimate the terms in the right hand side above. 

2. Estimates for the terms i?i and R^. 

The first and the last term are readily estimated. By Schwarz inequality and the 
definition (2.1) of the norm, we have 


(2.14) Ri = (V(u - u),\Iw)^ < ||u - u||/i||ic||?i. 

Schwarz inequality, the discrete trace inequality (2.2), and Lemma 2.4 give 


(2,15) 


EeGh ^ EeGh 


Eeg 
< Cdlrt- 


91 7 . -U nQPi f I 4- I 


E^Qh 


3. Splitting the boundary. 
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To treat the two remaining terms R 2 and i?3, we have to separate the cases when 
the edge size He is smaller or greater than e. To this end we denote the collection 
of edges of size greater than e by 

(2.16) Ql = {E&gh\e< /lu}, 
and the corresponding part of the boundary by 

(2.17) r, = IJ £;. 

E&G‘h 

We then write 

dv 1 1 

R2 + R3 = -{^--{uo-u)-g,w)^+ g + 

E&Gh 

(2.18) = ^ + 

EeGh ' 

= E {-+ 

EgTc 

EcFXFe 


4. Estimation of the contribution to R 2 + R 3 from the part Tg. 

On i? C Tg it holds e < hE and we estimate as follows, using Lemma 2.4, 

- E ^ ll^r(i')||Q,£ - (e + 


Ecr, 


(2.19) 


< X! I ^ H77r(T)||o,E • (e + hg)~^^^||w||o,E 


s(E 


2 hE 




(e + 7 /ib)2 


1/2 


EcTe ' ' EcTe 

< C{\\u - u|U + osc(/) + osc(e, uq, g)) || w|U. 


Next, we have to estimate 


( 2 . 20 ) 

We substitute 

( 2 . 21 ) 

which gives 

( 2 . 22 ) 


, dv 


-{uq -u)- g,w)^^. 


-{uo -u)+g 


du 


, dv 

^dn 


-{uq - u) - g,w)^^ = ( 


du 

dn 


dv 

dn 
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Now we define the strip 

(2.23) a = U K. 

KgTh 

Knr,^$ 

Following [1, 2, 7] we construct a linear finite element extension £hW G 14 of 
such that 

(2.24) £hw\r,=w\r, 
and 

(2.25) £hW = 0 in 
In [1, 2, 7] the following estimate is derived 

(2.26) llVarcIloA <C{J2 Ikllo.^) 

Ecr, 

We denote 

(2.27) r+ = anr. 

and split the boundary of flj in three parts (cf. Figure 2) 

(2.28) 9 a = r,u{r+\rju{ 5 a\r+}. 

Note that E^w ^ w on F+ \ F^. Since £hw\QQ^ \^p+ = 0, scaling and the estimate 
(2.26) show that 

(2.29) ( E + ( E 

E(zQ.^\rt 

< C\\V£hw\Wn, < C( E 

Ecr, 

and also 

(2.30) ( E ^E'll^/.^llo,E)'^Ec'||Vaic||oA <C( E ^e'I^IIE)'^'- 

Ecr+\r^ Ecr^ 

Further, since e < it holds 

(2.31) ( E < ^( E 

17'/—•n IT/—T' ^ 


\\w\\IeY^^<V2\\wU. 
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Next, integrating by parts and using (2.29)-(2.31) we estimate as follows 


(2.32) 


,du dv . 

= X! [-{f + + {V{u-v),V£hw) 


jfcn. 




2 \l /2 


tdn 

Ecn,\rt 

E C'( ^2 ^R'll/+ ^^IIo,a) ^ ( ^2 

Kcflc Kcflc 

+ ||V(m - u)||o,nJ|V£’ft,w||o,n, 

E<zn^\rt E(zn^\Tt 

-^{( ^aII/+ ^^IIo.r) ^ + llV(w — u)||o,n. 


Ken, 




Een,\Tt 


r9i;i 

-dn. 




W\\h. 


From a posteriori error analysis [3, 9] we know that 


(2.33) 


( Y 

Ecn,\r+ 


-dn. 



E C'( Y^ ^i'll/+ ^^IIo,a) ^ 

Ken, 


l|V(ri-n)||oA 


and 

(2.34) ( ^ hUf + ^vWIkY^ < C{\\\7iu - i;)||o,a + osc(/)). 

Ken, 


Hence we have 


(2.35) - ^,Ehw)^+ < c(^\\u-v\\h +osc(/))||r;||?,. 


Since £hW = rc on Fg, we get 

" 1^’ ^ C'(||u - vU + osc(/)) lliclU - (|^ - 
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For E C r+ \ Fe it holds that He < e < CHe with a constant that only depends on 
the regularity constant k. Thus we can estimate 

,du dv > 


dn dn 

= --{Rr{v),£hw)E - {u-v,£hw) 


E 


-Ecr+\r, 


(2.37) <C( ^ 


Bcr+\r, 


IE 


-||i?r(i^)||o,E||£/iR’||o,E 




<c( i: 


He 


Ecr+\r, 


{hs + tY 




0,E 




He + e 


E 




Ecr+\r, 


He + e 


Thus we have 


(2.38) (— - ^ C(||u- ■ullft, +osc(/) +osc(e,uo,5))lklU, 

which together with (2.19) gives 


s {-<1; - )<“» - 


1 


Ecr, 


'dn e ^’“/E ' 

(2.39) < C{\\u - + osc(/) + osc(e,Mo,5))||w|U. 

5. Estimation of the contribution to R 2 + R 3 from the part F \ Fg. 
It now holds e>hE- First write 

(2.40) 


dv I dv 1 1 

-- {uo-u}- g = - - [uo - V) - g + -(u - V) 

dn e dn e e 


= -Rr{v) + -{u — v). 


Hence, on E it holds 


- “> -S'”)* + 

--){Rr{v),w)^--{u-v,w)^ 


(2.41) 


< 


^e + y/iB e 
(g + ^/^^)gll-^r(^^)||o.E||w||o.E + -||m - 'c||o.e||R'||o.E- 


Since e + /i_e < 2e, it holds 


1 2 

-||m - 'c||o.e||w||o.e < , , ||m - wIIo.eIIwIIo.e- 

e e + hE 


(2.42) 













ERROR ANALYSIS OF NITSCHE’S METHOD FOR ROBIN BOUNDARY CONDITIONS 9 


Since /i£;/e < 1 we estimate as follows 


(2.43) 


ihE 

(e + 7/i£;)e 


||-Rr(i^)||o.£;||w||o.£; 


7/t^(e + 

(e + 7li£;)e 


1/2 


1/2 


= 7- 


T/2 (•£_!_ 


h 


1/2 


h 


1/2 


= 7- 


:l/2 


■7^e) 


l|/?r(^i)||o.E • (e + ^^^||w||o,e 

-—-ll/?r(^’)||o,E • (e + hE)~^^'^\\w\\o,E 

(1 + -^)^^^||i?r(ii)||o,E • (e + 


- + 7/tg) ^^^Ilii’llo.E- 

Combining (2.41)-(2.43) yields 

(2.44) ^ (^_{^^_\uo-u)-g,w)^ + j^^{Rr{v),w)^) 

Ecr\r, ' 

< C{\\u - v|U + osc(/) + osc(e, uq, g)) || w|U. 

6 . Collecting the estimates. 

Adding (2.39) and (2.44) gives 

(2.45) /?2 + /?3 < C'diu - ii|U +osc(/) + osc(e,Mo,5))lk||/i. 

The assertion then follows from this and (2.13), (2.14), and (2.15). 


□ 
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